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Abstract 


In studying the Earth’s geomagnetism, it has always been a challenge to separate magnetic fields 
from external currents originating from the ionosphere and magnetosphere. While the internal 
magnetic field changes very slowly in time scales of years and more, the ionospheric and 
magnetospheric current systems driven by the solar wind -magnetosphere interaction are very 
dynamic. They are intimately controlled by the ionospheric electrodynamics and ionosphere- 
magnetosphere coupling. Single spacecraft observations are not able to separate their spatial and 
temporal variations, and thus to accurately describe their configurations. To characterize and 
understand the external currents, satellite observations require both good spatial and temporal 
resolutions. This paper reviews our observations of the external currents from two recent LEO 
satellite missions: Space Technology 5 (ST-5), NASA’s first three- satellite constellation mission 
in LEO polar orbit, and Communications/Navigation Outage Forecasting System (C/NOFS), an 
equatorial satellite developed by US Air Force Research Laboratory. We present 
recommendations for future geomagnetism missions based on these observations. 
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1. Introduction 


The existence of the Earth’s internal magnetic field is vital to life on Earth because it acts 
as a giant shield to protect the Earth from the solar wind (charged particles from the sun) and 
cosmic rays. The Earth’s main magnetic field is generated by an internal electric current 
maintained by a rotating and electrically conducting fluid in the Earth’s outer core powered by 
the convective geodynamo [Glatzmaier and Roberts, 1995a, b]. Studying the geomagnetism by 
measuring and monitoring the Earth’s magnetic field provides an important way to probe the 
Earth’s liquid core and its change with time. However, other sources of magnetisms, although 
small in comparison with the main field from the internal source, also contribute to the Earth’s 
magnetic field. They include crustal magnetic fields, ocean currents, and external currents 
originating from the ionosphere and the magnetosphere. At any location and any moment, the 
magnetic field is the vector sum of the fields from all these sources. External currents originating 
from the ionosphere and the magnetosphere contaminate the magnetic field measurements in 
geomagnetism. In studying the Earth’s geomagnetism, it has always been a big challenge to 
separate the magnetic fields from internal and external sources. Although there is a long history 
of direct measurements of the Earth’s magnetic field on the ground, ground-based observations 
suffer from uneven spatial coverage and large gaps in the oceans. Only Low Earth Orbit (LEO) 
satellites can provide true global mapping of the Earth’s magnetic field. High precision 
measurements from dedicated geomagnetism satellites such as Magsat (1979-1980), Orsted 
(since 1999), and CHAMP (2000-2010) have resulted in significant advances in monitoring, 
modeling, and understanding the Earth’s magnetic field [see Olsen and Stolle, 2012 for a 
review]. 
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Geomagnetic field changes on various time scales. The internal magnetic field changes 
very slowly, in time scales of years and more. But the external currents are very dynamic and 
vary in much shorter time scales (seconds to days). To characterize and understand the external 
currents, satellite observations require both good spatial and temporal resolutions. Single 
spacecraft measurements do not allow us to separate spatial and temporal variations, and thus are 
unable to accurately and fully describe their configurations. This paper will discuss our recent 
LEO spacecraft observations of the external currents driven by the solar wind-magnetosphere 
interaction and recommendations for future geomagnetism missions. 
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2. External Current System Driven by the Solar Wind-Magnetosphere Interaction 

To a first order approximation, the Earth’s internal magnetic field is dipolar resembling 
that of a bar magnet. It acts as an effective obstacle to the flow of charged particles from the sun, 
called the solar wind. The interaction with the solar wind flow confines the Earth’s magnetic 
field in a cavity called the magnetosphere by compressing the dayside and stretching the 
nightside magnetic field lines. The magnetospheric cavity has a compressed dayside and a long 
comet-like tail, which is significantly distorted from the dipolar magnetic field. Such a distortion 
and the overall shape of the magnetosphere are the direct result of the presence of large-scale 
electric currents systems in the magnetosphere and the ionosphere that are driven by the solar 
wind-magnetosphere interaction. Figure 1 illustrates the large-scale electric current systems in 
the magnetosphere and the ionosphere. They include: the magnetopause (Chapman-Ferraro) 
current flowing on the magnetosphere boundary, the ring current in the inner magnetosphere, the 
tail current flowing in the neutral sheet across the magnetotail, field-aligned (Birkeland) currents 
flowing in and out of the ionosphere and coupling the magnetosphere to the ionosphere, as well 
as associated horizontal currents in the ionosphere. The horizontal ionospheric currents include 
Pedersen currents in the auroral zone and across the polar cap, and auroral electrojets (Hall 
currents) around the auroral oval. These magnetospheric and ionospheric currents respond 
dynamically to variations of the solar wind plasma and the interplanetary magnetic field. 
Changes of these current systems cause geomagnetic disturbances. Thus, the solar wind- 
magnetosphere interaction is the main driver for various geomagnetic activities in short time 
scales (seconds to days), and the ultimate energy source is provided by the solar wind from the 
Sun’s atmosphere. The term “space weather” is used to refer to the changing environment of 
plasma, magnetic fields, and radiation in near-Earth and interplanetary space due to solar 
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variability. Large-scale electric currents in the magnetosphere and ionosphere constitute 
important space weather parameters. During magnetic storms and substorms, these currents 
intensify in response to the enhanced solar wind-magnetosphere interaction. 

In the region below ~ 1000 km from the surface of the Earth, where geomagnetism 
satellites fly, the external currents that generate the largest magnetic fields in contaminating the 
geomagnetism measurements are field-aligned currents (FACs) at auroral latitudes, horizontal 
currents in high latitude ionosphere, and the ring current in the inner magnetosphere. Among 
them, field-aligned currents flow into and out of the ionosphere in the auroral zone, and are 
closed by horizontal Pedersen currents to complete the current loops in the auroral zone and 
across the polar cap in the ionosphere. The combined FAC-Pedersen current loops are mostly 
invisible on the ground because the magnetic fields are confined within the current loops. But 
polar orbiting LEO satellites pass right through the FAC layers and make direct in situ 
measurements of the magnetic field disturbances generated by the combined FAC-Pedersen 
current loops. Their magnetic field perturbations are transverse to the background magnetic field, 
and can reach to over 1000 nT in the magnetic field components. There is no perturbation in the 
magnetic field strength in the in-situ measurements because the effect of the transverse magnetic 
field perturbation is to twist the magnetic field lines without changing the field strength. Auroral 
electrojects (Hall currents) are another type of large-scale horizontal currents in high latitude 
ionosphere. They flow in the auroral oval, westward in the dawn side and eastward in the 
duskside. They are largely closed within themselves in the polar ionosphere. During substorms, 
enhanced westward auroral electrojets, called the substorm current wedge, are fed by field- 
aligned currents from disrupted tail current [McPherron et al., 1973], The magnetic field 
signatures of auroral electrojets can be readily measured on the ground below the auroral zone 
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[e.g., Kamide et al., 1981; Friis-Christersen et al., 1985], In space, auroral electrojets flow below 
LEO satellites. They cause a magnetic perturbation mainly in the magnetic field strength in 
auroral latitudes, either positive or negative depending on the local time [Zanetti et al., 1984], 
which decreases with the altitude and can be detected at altitudes below -700 km [Moretto et al., 
2002; Le et al., 2009], 

In the magnetosphere, the ring current, the tail current, and the magnetopause current are 
all remote current systems to LEO satellites; and they all produce global magnetic disturbances 
that can be readily measured on the ground and by LEO satellites. However, the magnetic fields 
associated with the magnetopause current and the tail current are relatively small because they 
are more remote. The ring current in the inner magnetosphere makes the most significant 
contribution to these global disturbances due to its relative proximity to the Earth. The ring 
current is formed by charged particles in the magnetosphere that are trapped in the Earth’s 
magnetic field from the solar wind through enhanced solar wind-magnetosphere interaction. It 
flows westward in the equatorial magnetosphere and produces a global southward magnetic field 
perturbation at the Earth. Since the Earth’s main magnetic field is northward in the equatorial 
region, the ring current causes a global depression of the magnetic field strength, and the 
equatorial average value of which, the Dst index, is used to monitor and characterize the ring 
current. It has been showed that the absolute value of the Dst index is proportional to the total 
energy content of the charged particles in the ring current region [Dessler and Parker, 1959; 
Sckopke, 1966], A prolonged negative Dst index is an indication of a magnetic storm in progress 
and also a measure of the storm intensity. The more negative the Dst index is, the more intense 
the magnetic storm is. Geomagnetic storms are classified based on the Dst index, as moderate 
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(Dst > -100 nT), intense (-250 nT < Dst < -100 nT) and super-storm (Dst < -250 nT). During the 
solar cycle 23 (May, 1996 - November, 2008), 11 super-storms occurred [Echer et al., 2008], 

Both the magnetopause current and the tail current add to the magnetic disturbances at 
the Earth and contribute to the Dst index, although to a much lesser extent. The magnetopause 
current is the boundary of the Earth’s magnetic field. It flows from dawn to dusk in the dayside, 
opposite to the ring current. It is controlled by the solar wind dynamic pressure, and contributes a 
northward magnetic field of ~ +20 nT at the Earth and an positive value of the same amount in 
the Dst index during nominal solar wind conditions [Burton et al., 1975], The Dst index also 
shows a +20-30 nT sudden rise, the so-called the storm sudden commencement (SCC), in 
response to a sudden increases of the solar wind dynamic pressure at the beginning of a classic 
magnetic storm [Dessler et al., I960], The tail current, on the other hand, flows from dawn to 
dusk in the night side in the same sense of the ring current. Its contribution to the Dst index can 
be significant during storms and substorms. Observations show that the tail current can account 
for ~ 20-25% of the measured Dst index variation during storms and substorms [Turner et al., 
2000; Ohtani et al., 2001], 

Although the Dst index is an indication of the ring current strength, it does not provide 
any information about the local time asymmetry of the ring current. Both ground-based and in 
situ satellite observations have provided evidence that the ring current has both a symmetric and 
asymmetric parts, especially during stormtime [Fukushima and Kamide, 1973; Iyemori, 2000; 
Greenspan and Hamilton, 2000; Turner et al., 2001], A significant fraction of the ring current is 
partial, which flows only within a limited longitudinal region and must be diverted out of the 
equatorial region as field-aligned currents to close in the ionosphere. The ring current 
distributions deduced from in situ magnetic field data show that the partial current is much 
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stronger than the symmetric current, up to a factor of 5 under moderate storm conditions [Le et 
al., 2004], Thus the partial ring current makes the major contribution to the Dst index. To 
describe the asymmetric nature of the ring current, a new set of geomagnetic disturbance indices, 
a longitudinal asymmetric (ASY) and a symmetric (SYM) indices, are introduced for both the H 
and D components of the magnetic field at the surface of the Earth at mid-latitude [Iyemori et al., 
1992], The SYM-H is essentially the same as the hourly Dst index but with higher resolution, 
which is the average disturbance at every minute for the H-components of all the stations. The 
partial ring current contributes to the SYM-H index the same way as it does to the Dst index. 
Thus, SYM-H does not represent the strength of the symmetric ring current, but the average ring 
current strength for both the symmetric and asymmetric components. On the other hand, the 
ASY-H index, which is the range between the maximum and minimum deviation of the H- 
components from the SYM-H, is an indication of how asymmetric the ring current is. 

In summary, external currents in the magnetosphere and ionosphere are very dynamic 
and respond to the solar wind-magnetosphere interaction. Separating these external currents from 
geomagnetism measurements requires characterizing their strength, spatial variation, and 
temporal evolution for both quiet and disturbed times. The main goals of external current 
investigations are to understand how they vary with solar wind parameters, how they vary with 
location and local time, and how they change with time. While a single polar-orbiting LEO 
satellite covers all latitudes for two local times once every ~ 90 minutes, it cannot separate 
spatial and temporal variations. It is desirable to have significantly denser coverage in space and 
time with a multi-satellite constellation. The upcoming ESA’s magnetic field mission Swarm 
will be the first constellation of satellites for geomagnetism and is expected to lead to new 
insight into many natural processes responsible for the Earth’s magnetic field, including the solar 
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wind-magnetosphere interaction [Friis-Christensen et al., 2006], In the following section, we will 
review our recent observations of external currents from two LEO satellite missions. One of the 
missions is Space Technology 5 (ST-5), NASA’s first three-satellite constellation mission in 
LEO polar orbit as shown in Figure 2 [Slavin et al., 2008], The other one is 
Communications/Navigation Outage Forecasting System (C/NOFS), an equatorial satellite 
developed by US Air Force Research Laboratory [de La Beaujardiere, 2004, 2009]. Although 
these satellites were not equipped with instruments for geomagnetism purposes, they all carried 
research-quality magnetometers for studying the external currents in the ionosphere and the 
magnetosphere. These two missions have provided us important magnetic field data for 
understanding time-space characteristics of the external currents and valuable lessons for 
designing post-Swarm geomagnetism missions. 
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3. ST-5 Observations of Field-Aligned Currents and Ionospheric Currents 

ST-5 is a three micro-satellite constellation deployed into an elliptical (300 km perigee 
and 4500 km apogee), dawn-dusk, sun-synchronous polar orbit from March 22 to June 21, 2006, 
for technology validations. The three spacecraft are maintained in a string-of-pearl constellation 
with controlled spacing ranging from under 50 km up to ~ 5000 [ref. Figure 1 in Slavin et al., 
2008] . Each spacecraft carried a boom-mounted miniature tri-axial fluxgate magnetometer, and 
returned high quality magnetic field data as the constellation flew in formation and made 
simultaneous multi-point measurements of the magnetic field through the Earth’s dynamic 
ionospheric current systems. A substantial volume of magnetic field data was taken over a range 
of inter-satellite spacing. These separations allow us to determine the properties of FACs and 
separate spatial versus temporal structures of auroral field-aligned currents over a wide range of 
spatial (~ 50-4000 km) and temporal (~ 5 s-10 min) scales. 

Field-aligned currents usually appear as quasi-planar “sheets” that tend to be loosely 
parallel to lines of constant geomagnetic latitude [e.g., Iijima and Potemra, 1978], Typically, 
there is a set of “Region 1” or “Rl” FACs along the high latitude edge of the auroral oval, which 
originate near the equatorial edge of the magnetosphere. The Rl currents flow into the 
ionosphere in the dawnside and out of the ionosphere in the duskside. At the lower latitude edge 
of the auroral oval there is also a set of “Region 2” or “R2” FACs with polarities opposite to Rl , 
which originate in the region where the ring current has a divergence due to the existence of a 
partial ring current. The interaction between the solar wind and the magnetosphere is controlled 
by the interplanetary magnetic field (IMF) and solar wind conditions [e.g., Cowley, 1984], The 
IMF and solar wind constantly change, making the field- aligned current systems highly dynamic. 
Temporal variability of the field-aligned currents at time scales less than the orbit period of low 
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Earth orbit spacecraft (~ 90 min) cannot be assessed using data from single spacecraft. The data 
from ISEE 1 and 2 magnetometers provided the first dual-point simultaneous measurements of 
FACs at mid-altitudes (2.4-7 R E ). The four- spacecraft Cluster data have also been used to study 
FACs at mid- and high altitudes (4-11 R E ) [e.g., Cargill et al., 2001; Johansson et al., 2004; 
Draper et al., 2005; Figueiredo et al., 2006], The 3-spacecraft ST5 mission provides the first 
multi-point measurements of FACs at low altitudes (~ 300 - 4500 km), which are 
complementary to the mid- and high-altitude observations. 

3.1. The current density, motion and velocity of FACs 

Previously, the standard method for calculating the current density from single spacecraft 
magnetic field data requires the assumption that the FAC is a stationary, infinite current sheet to 
the east-west direction [Iijima and Potemra, 1976], As the spacecraft passing through the 
stationary current sheet, the magnetic field perturbation in the eastward component 6B E and the 
spacecraft velocity in the northward direction V s/C (same as the current sheet normal) are used to 
calculated the current density: 

J n ={-llli 0 )(\IV slc ){dB E ldt) (1) 

Similarly, the thickness of the current sheet can be determined as L= V s/c -dt, where 8t is the time 
duration of the current sheet crossing. The errors in such calculations are positively correlated to 
the ratio of W cs /V s/c \, where V cs is the northward velocity component of the current sheet motion. 
The simultaneous multi-point measurements from ST-5 constellation allow us to determine the 
velocity of the current sheet motion and thus to relax the assumption that the current sheet be 
stationary. From the locations of the spacecraft and the times when two spacecraft observe the 
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same current sheet structure, we can determine V cs . Then the current density can be determined 
more accurately by correcting for the motion of the current sheet [Slavin et al., 2008]: 

j„ =(-i/^)[i/(v s/c -v cs mdB E /dt) (2) 

The statistical study using the entire ST-5 data set shows that the current sheet velocity is 
quite variable and occurs in a large range from -1 to 1 km/s at ST-5 altitudes of ~ 300-5000 km; 
and current sheets tend to move faster/slower during intervals of higher/lower geomagnetic 
activities [Wang et al., 2009]. The ratio \V CS /V S/C \ occurs in the range ~ 0 - 25% with the median 
(mean) value of 4% (6%). The large range of W cs /V S i C \ happens for all time periods with both 
high (Kp >4) and low (Kp < 4) geomagnetic activities. During periods of low geomagnetic 
activities, there is still a significant fraction of the events with the ratio I V cs /V s/c \ higher than 
10 %. 


The ST-5 multi-point measurements of FACs also allow us to measure the current density 
using the gradiometry technique pioneered by the 4-spacecraft Cluster mission [Balogh et al., 
1997], It is the first mission to provide the necessary multi-point measurements to support 
magnetic gradiometry in low Earth orbit [Slavin et al., 2008], When two spacecraft are within a 
current sheet simultaneously, the current density can be determined by the gradient of the 
magnetic field measured at the two spacecraft. This method has the advantage of removing 
contaminations due temporal variations in the calculation. Temporal variations with wavelengths 
comparable to or greater than the spacecraft separation (e.g., Alfven waves) are measured 
simultaneously by the two spacecraft and thus removed in computing the gradients. ST-5 
provided numerous opportunities for applying the gradiometry technique when the inter- 
spacecraft separations went down to ~ 100 km or less. 
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Figure 3 is adapted from Slavin et al. [2008] showing two examples of FAC current 
density determination by the magnetic grodiometry for auroral oval passes on June 15 and 20, 
2006. In ST-5 constellation, the leading, middle, and trailing spacecraft are named 155, 094, and 
224, respectively, in the string-of-pearl configuration. For the June 15, 2006 auroral oval 
crossing the 094-155 and 094-224 s/c separations were 131 and 100 km, respectively. The 
thicknesses of the R1 and R2 intervals were 524 and 442 km, respectively. Similarly, for the 20 
June 2006 event, the 094 - 155 and 094 - 224 separations were 49 and 92 km, respectively. The 
thicknesses of the R1 and R2 intervals were 176 and 223 km. Thus in both cases, the spacecraft 
separations are well below the current sheet thickness; the data from both the leading pair 094- 
155 and the trailing pair 094-224 can be used to calculate the magnetic gradient. Figure 3 
displays the FAC current density using the gradiometry technique for the leading and trailing 
pairs along with the SC094 single spacecraft, motion-corrected current density. For the June 15, 
2006 FAC event (left panel of Figure 3), all three traces are very similar. Hence, we conclude 
that for this interval the FACs were very stable over a time span of several minutes. On the other 
hand, there are differences of up to a factor of 2 between the peak current densities determined 
from leading and trailing pairs for the June 20, 2006 FAC event. The single s/c current density 
determination from SC094 only agrees well with the values determined by the trailing pair. Since 
the leading spacecraft 155 was only ~ 10 s ahead of the middle spacecraft 094, there are large, up 
to 50% changes took place during this brief interval. 

3.2. Spatial and temporal variability. o£FACs 

Field-aligned currents not only are in motion, but also change with time. Single spacecraft 
measurements are unable separate their spatial and temporal variations. Temporal variability in 
time scales less than ~ 100 min (the orbit period of low Earth orbit spacecraft) cannot be 
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assessed using data from single spacecraft. The data from ST-5 constellation provide the first in 
situ observations of FAC temporal variability at low altitudes in time scales of -10 min and less. 
As the three spacecraft cross the FAC region successively along the same trajectory, their 
magnetic field profiles would exactly track each other with only time delays when the magnetic 
variations are due to spatial changes. But any differences in the magnetic field profiles would 
indicate temporal changes of the current sheet structures. Thus, we can study the temporal 
variability of the FACs using the magnetic field profiles from multiple spacecraft in a string-of- 
pearl configuration. 

Figure 4 is adapted from Le et al. [2009] showing an overview of the northern and 
southern polar cap passes during an intense magnetic storm on 14 April 2006. The Dst index at 
its minimum is -111 nT. The spacecraft trajectories are shown in Figure 4 (top), in which we 
have mapped the spacecraft positions to their ionospheric footprints at 110 km altitude along the 
magnetic field lines. The nested circles represent constant magnetic latitudes separated by 10° 
and centered at the Earth’s magnetic Pole. (Note we have flipped the trajectory for the southern 
polar cap pass to the northern hemisphere.) The trajectories of the three ST5 spacecraft 094, 155 
and 224 are color coded in black, red and blue, respectively. The tick marks on the trajectories 
are separated by 10 min with color-coded time tags. During this time period, the spacecraft move 
across the polar cap near the dawn-dusk meridian plane: from dusk (dawn) to dawn (dusk) in the 
northern (southern) hemisphere. In the string-of-pearls configuration, the middle spacecraft 094 
and the trailing spacecraft 224 are close together and have a large separation from the leading 
spacecraft 155. The lag time along the orbit is about 10 min 094-155 pair and about 1 min for 
224-094 pair. The spacecraft spacing is -5331 (-3564) km for 155-094 pair and -531 (-368) km 
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for 094-224 pair over the northern (southern) polar cap. Thus, the observations from this pass 
allow us to evaluate the FAC variability at these two temporal and spatial scales. 

The bottom panels of Figure 4 shows an overview of ST5 magnetic field variations 
generated by FACs during these two passes, including the three components of the magnetic 
field residual vector (data with the internal IGRF model magnetic field removed) in the solar 
magnetic (SM) coordinate system, as well as the residual of the magnetic field strength. The data 
from the three spacecraft are also color-coded, but the labels for the spacecraft positions 
(altitudes, magnetic latitudes and magnetic local times) are for the middle spacecraft 094 only. 
Since it is an active period, we have observed strong FAC activities in the auroral region, in both 
dawn and dusk sides, as evident by the perturbations of magnetic field components as large as ~ 
1000 nT in the bottom panels of Figure 4. Since FACs flow along the background magnetic field 
direction, their magnetic field perturbations are transverse to the background field with the 
magnitude much smaller than the background field strength. The presence of FACs only twists 
the field line direction without changing the field strength, and thus, we do not expect any 
perturbations in the magnetic field strength. However, we do see some perturbations of the 
magnetic field strength associated with FACs in the northern polar cap passes near the orbit 
perigee with altitudes less than ~ 1000 km, such as the case in the bottom left panel of Figure 4. 
These perturbations in the field strength are generated by strong auroral electrojets in the auroral 
zone during this active interval. As the altitude increases, their magnetic perturbations decrease 
and eventually become undetectable, as in the case in the bottom right at altitudes ~ 3500 km. 

In order to examine the temporal variability of FAC structures, Figure 5 displays the time 
shifted magnetic field data from the three spacecraft for the April 14, 2006, northern polar cap 
pass near the perigee, which is in the altitude range of interest to geomagnetism missions. It is 
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also adapted from Le et al. [2009]. In this figure, the magnetic field data from the leading (red 
for SC 155) and trailing (blue for SC224) spacecraft are time- shifted to match the middle 
spacecraft (black for SC094) in order to line up the large-scale current structures observed by the 
three spacecraft. The time shifts are determined from the cross-correlation analysis of large-scale 
current structures in 094-155 and 094-224 pairs. The time labels on the bottom of the horizontal 
axis are for SC094. The magnetic field components are displayed in the coordinate system 
determined by the Minimum Variance Analysis (MVA) [Sonnerup and Cahill, 1967], where i,j, 
and k are the maximum, intermediate, and minimum variance directions, respectively. For FACs, 
k is also the field-aligned direction, and i and j are the two directions transverse to the magnetic 
field. If the magnetic field variations are mainly seen in one component (8Bi), the infinite current 
sheet would be a good approximation for FACs, and the current sheet normal direction would be 
along the k direction. Otherwise, the magnetic field variations would be seen in both i and j 
directions for a current sheet with finite width. 

For the duskside FACs (left panel of Figure 5), large-scale magnetic field variations are 
mainly in the &B t component, and the bB i component contains mainly mesoscale variations with 
much smaller amplitudes at all three spacecraft. These observations indicate that the infinite 
current sheet approximation applies to the large-scale FACs observed here. Since the mesoscale 
variations, which are embedded within the large-scale currents, are seen in both t>B [ and 82?j, the 
mesoscale currents are filamentary, and cannot be treated as infinite sheet currents. For the 
dawnside FACs in the right panel, the large-scale FACs still appear mainly in the h>B t component, 
but there are significant mesoscale variations in the &B } component at all three spacecraft. The 
amplitudes of the h>B } variations become comparable to those in t>B [ in this case. Clearly the 
mesoscale currents are generally in the form of current filaments. The traditional method for 
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determining the current density described in Equations (1) or (2), which has assumed the infinite 
current sheet approximation, will not be applicable to the mesoscale FACs. In this case, 
calculating the current density requires the knowledge of both bB, gradient in the j direction and 
bB } gradient in the i direction. Although the three ST-5 spacecraft in string-of-pearl enable us to 
study the temporal variability of the currents, they are not in the most desirable configuration for 
measuring the density of these mesoscale currents because they do not provide adequate 
separations in the east- west direction (the i direction). It is most desirable that three spacecraft 
are in a triangular configuration in the plane perpendicular to the magnetic field. 

From the time- shifted magnetic field data from the three spacecraft in Figure 5, the 
characteristics of FACs temporal variations are very evident for the both FAC intervals. First of 
all, large-scale currents, shown mainly in bB r are relatively stable in time scales of ~ 10 min. 
The duration and magnitude of the overall bB t variations maintain the similar values at the three 
spacecraft. However, the mesoscale structures embedded within the large-scale currents show 
significant changes in the same time scales. Comparing the observations of both bB t and bB } 
components of ~10 min apart (red versus black/blue traces), the mesoscale structures at SC155 
exhibit the largest differences from those of SC094 and SC224. We can observe changes in 
magnitude, polarity, as well as locations for the mesoscale currents. Meanwhile, the data also 
show that the time scales for the currents to be relatively stable are ~1 min for mesoscale 
currents and at least ~ 10 min for large-scale current sheets. 

3.3. Ionospheric closure of FACs 

Pedersen currents in the ionosphere are the closure currents for FACs. The combined 
FAC-Pedersen current loops are shown in Figure 1. Near the dawnside (dusk) auroral oval. 
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region 1 FACs flow into (out of) the ionosphere at the high-latitude edge of the oval; region 2 
FACs flow out of (into) the ionosphere. Most of the current closure takes place via local 
Pedersen currents within the auroral zone flowing between the upward and downward FAC pair, 
i.e., Pederson currents flow equatorward (poleward) at the dawnside (duskside) auroral zone to 
form a closed current loop. However, observations show that there is generally an imbalance 
between the R1-R2 pair in either dawnside or duskside, i.e., the total current flowing in R1 is 
more than that in R2 [Iijima and Potemra, 1976; Sugiura and Potemra, 1976], Thus, there are net 
currents into (out of) the ionosphere due to the R1-R2 imbalance in the dawnside (duskside) 
auroral region. Such net currents need to be closed within the R1 FACs on either side of the pole 
via cross-polar cap Pedersen currents, also shown in Figure 1. 

We can use a simplified model to calculate the magnetic field perturbations expected 
from the combined field-aligned current-Pedersen current system. Figure 6 is adapted from Le et 
al. [2010] showing the FAC current setup and geometry for simple calculations of the magnetic 
field signatures. The simplified geometry is such that the X direction is from dawn to dusk with 
the magnetic pole at X = 0, Z is vertically up along the magnetic pole, and Y points into the 
paper, westward (east- ward) in the dawnside (duskside). The infinite planar current sheets are in 
the YZ plane with current flowing directions shown as arrows in Figure 6a. The three pairs of 
balanced current sheets in Figure 6a (left) are equivalent to the two pairs of unbalanced current 
sheets in Figure 6a (right). In Figure 6b, we first calculate the magnetic field from two pairs of 
balanced R1-R2 currents on each side of the pole using characteristic current properties listed in 
the left panel. In this case, the R1 and R2 are balanced and the net current on either side of the 
magnetic pole is 0. The calculated magnetic field in Figure 6b (right) is the well-known unipolar 
bump in the azimuthal direction (the Y direction) on either side of the magnetic pole. The east- 
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west component of the magnetic field 6 B y is mainly confined within the R1-R2 current sheets 
and quickly decreases to zero away from the current pair, both over the pole and equatorward 
from the R1-R2 currents. Next, we decrease the current intensity of the R2 current by 25% so 
that the R1-R2 currents are imbalanced, as shown in Figure 6c (left). The net current flowing into 
(out of) the ionosphere is 25% of the total R1 current in the dawnside (duskside). The magnetic 
field bB y within the R1-R2 circuit remains to be unipolar with reduced magnitude as shown in 
Figure 6c (right). But there appears to be a magnetic field offset over the pole between the 
dawnside and duskside FACs. If we further decrease the R2 current intensity so that the net 
current is 50% of the total R1 current, the magnetic field bB y offset over the polar cap also 
increases, as shown in Figure 6d. Thus, the signature of the imbalanced R1-R2 pairs is the 
magnetic field offset over the polar cap. Although the actual FACs and ionospheric current 
systems are much more complex than this simple model illustrates, it demonstrates the type of 
magnetic signatures and their magnitudes we expect to observe in situ. Using this offset, we can 
quantify the R1-R2 imbalance based on in situ magnetic field observations from polar-orbiting 
spacecraft. 

Figure 7 displays two examples of ST-5 polar cap crossings showing the magnetic field 
observations and the deduced current density along the orbit track (adapted from Le et al. 
[2010]). The horizontal axis in each panel is the spacecraft distance from the magnetic pole. 
From both the examples, it is very clear that there are indeed magnetic field offsets in bB y across 
the polar cap (bottom panel), indicating that the R1 currents are stronger than the R2 currents and 
there are net currents flowing into or out of the ionosphere. In order to quantify the imbalance of 
the R1-R2 FACs, we calculate the total current intensity using the magnetic field observations 
for each pass and determine the net current (top panel). Ideally we would like to have two 
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spacecraft, one on either side of the pole, to measure the dawnside and duskside FACs 
simultaneously. Since the largest time lag of the three ST-5 spacecraft is only ~ 10 min, we do 
not have the cases when the dawnside and duskside currents are observed simultaneously. Thus 
we measure the net current density at the dawnside and duskside individually and examine them 
statistically. Figure 8 shows the scatter plot of the R2 current intensity versus the R1 current 
intensity in the duskside and dawnside, respectively, for all the ST-5 events. In each panel, the 
solid line has a slope that is the average of the R2 intensity to R1 intensity ratio. The dashed line 
has a slope of 1, where the R1 and R2 currents have the same intensity. In both the dawnside and 
the duskside, almost all the data points are located in one side of the dashed line, where the R1 
currents are stronger than the R2 currents. The net currents, due to this R1-R2 imbalance, are 
about 5% of the R1 currents on average in both sides of the pole. This net current will flow as 
Pedersen current across the polar cap in order to close the imbalanced FACs in the ionosphere. 
Although the cross-polar cap Pedersen currents are only a small fraction of the R1 currents, they 
still represent a significant amount of Pedersen currents flowing across the polar cap. Previous 
observations have determined that the total R1 currents are in the order of a few MA, comparable 
to the total amount of Chapman-Ferraro current in the magnetopause [e.g., Midgley and Davis, 
1963] and the ring current in the inner magnetosphere [e.g., Le et al., 2004], Thus, the total 
amount of the cross-polar cap Pedersen currents is in the order of ~0.1 MA. Despite the fact that 
the R1-R2 imbalance only contributes ~5% of the total R1 currents to the cross-polar cap 
Pedersen currents whereas ~95% flow as auroral zone Pedersen currents, the integrated Joule 
heating rate of the cross-polar cap Pedersen current accounts for a much larger fraction due to the 
much larger area they flow in the polar cap. Hence, the associated energy dissipation in the polar 
cap cannot be ignored. 
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4. C/NOFS Observations of the Ring Current During Magnetic Storms 


The C/NOFS spacecraft was launched into a nearly circular 13° inclined orbit on 17 April 
2008 with a scientific payload designed to specify and forecast plasma density irregularities in 
the equatorial ionosphere that degrade trans-ionospheric radio transmissions [de La Beaujardiere 
et al., 2004, 2009]. The single satellite is 3-axis stabilized and has an orbital period of ~ 97 min. 
Initial apogee and perigee were at altitudes of 867 and 401 km, respectively. The Vector Electric 
Field Instrument (VEFI) suite on the C/NOFS spacecraft includes a sensitive 3-axis fluxgate 
magnetometer mounted on a 0.6 m boom [Pfaff et al., 2010], Measurements yield full magnetic 
vectors every second over the range of ±45,000 nT with a one-bit resolution of 1.37 nT in each 
component. During magnetic storms, the ring current produces the dominant external magnetic 
field in the equatorial region. C/NOFS provides a complete coverage of all local times every ~ 
97 min, a time scale much smaller than the life span of magnetic storms. Thus, C/NOFS 
magnetic field measurements enable us to study local time variations of the ring current and its 
evolution during storms. Herein we demonstrate that a single equatorial LEO satellite enables us 
to monitor and track the ring current evolution, study the local time variation, and calculate near 
real time Dst index. 

In the low-latitude ionosphere, the ring current is expected to produce a negative 
perturbation in the northward magnetic component ( bB N ). Thus, we concentrate on bB N data 
observed during magnetic storms to examine the ring current characteristics. Figure 9 shows the 
IMF, the solar wind and the Dst index during the July 22, 2009 magnetic storm, which is one of 
the events studied in Le et al. [2011], The magnetic storm started shortly after the arrival of an 
interplanetary shock at ~ 01:00 UT on July 22. The main phase minimum of -79 nT in Dst was 
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reached at ~ 09:00 UT. This moderate magnetic storm was a consequence of the strong 
southward turning of IMF Bz after the shock compression. 

Figure 10 shows the local time variations of 5 Bn at the six stages of the storm’s 
development corresponding to the vertical dashed lines in Figure 9. Plots in Figure 10 use a 
format similar to that introduced by Love and Gannon [2010], In each panel, a full orbit of 5 Bn 
measurements are displayed as a function of the spacecraft magnetic local time (MLT). The 
baseline for 65 at= 0 is denoted by the dashed line circle. Positive/negative bB N is plotted 
inside/outside the baseline circle. Blue circles represent the Dst index as the radial separations 
between the blue and baseline circles is the Dst value. Solid black traces represent bB N plotted as 
a function of MLT. Similarly, the radial distance between black traces and baseline circles 
represent the values of bB N . Since the C/NOFS orbital plane does not align with the magnetic 
equator, bB N at the spacecraft is normalized by the cosine of the magnetic latitudes where 
measurements were made. Thus, displayed bB N is the component of the magnetic field residual 
parallel to the geomagnetic dipole axis. The red circle in each panel represents the least square fit 
to bBN using an off-center circle. The circle fitting results in two fitting parameters: the center 
and the radius of the fitting circle. Small red crosses mark the centers of fitted circles. The center 
of the fitted circle provides information about the local time asymmetry of bB N : (1) its MLT 
indicates where maximum bB N occurs; and (2) its radial displacement from the origin is a 
measure of the degree of the MLT asymmetry of bB N . Radii of fitted circles (after 100 nT 
baseline removal) are nearly identical to the absolute value of the orbital-averaged bB N - As an 
analogy to how the Dst index is estimated from the ground-based bB N , the orbital- averaged bB N 
from C/NOFS data can be used as a real-time provisional Dst index. We note that the equatorial 
electrojet (EEJ) also contributes to equatorial bB N in the dayside. On the ground, its signal 
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reaches up to 80 nT near the dayside magnetic equator [Manoj et al., 2006], and their effect is 
avoid by using mid-latitude ground stations in the calculation of the Dst index. Spacecraft 
observations show that EEJ signals are confined mainly within ±3° from the magnetic equator 
and maximize between 10:30-12:00 LT; and their bB N magnitudes are in the order of 20 nT at ~ 
450 km and ~ 10 nT at ~ 700 km [Aiken and Maus, 2007], Thus the EEJ magnetic signals near 
the magnetic equator are in the same order of the quiet time ring current in typical ionospheric 
satellite altitudes of ~ 400-700 km range, and much smaller than those of stormtime ring current. 
Here we ignored the EEJ effect in the study of the stormtime ring current as the spacecraft is 
outside the EEJ region for most of the orbit. 

Figure 10 provides a basis for perceiving the morphology of the ring current’s local time 
evolution. Panel 1 describes the pre-storm situation when the Dst and 8 Bn traces were very close 
to the baseline. Even in very quiet times, the red cross centroid was slightly shifted (4.1 nT) from 
the origin toward the pre-midnight sector, indicating the ring current is slightly asymmetric. 
Panels 2 and 3 show bB N distributions measured during the early main phase and at maximum 
epoch, when the centroid was shifted toward the dusk-evening MLT sector by 30.2 nT and 55.6 
nT, respectively. This is a sign that the storm time ring current quickly becomes very asymmetric 
during the main phase. Comparing the red and blue circles we see that near the dawn meridian 
(where the minimum bBN occurs) Dst was slightly more negative than bB,\i- However, similar to 
DMSP observations [Burke et al., 2011], at evening- midnight local times bBN was significantly 
more negative than Dst. This asymmetry is contributed by the rapid development of a partial ring 
current as well as the remote field-aligned currents that close the partial ring current in the 
ionosphere. The maximum of bB N is in the evening-midnight section during the main phase. 
Panel 4 is near the early recovery phase and a slight Dst dip to a second minimum the bB N 
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distribution appears to be far more symmetric than was detected during the two previous orbits. 
The displacement of the red cross centroid moved back to 10.9 nT. Thus, the ring current 
recovery in this case started with a rapid decay of the partial ring current. Panels 5 and 6 indicate 
that during the later parts of the recovery phase the Dst and 6 B N distribution traces come closer 
together suggesting that the ring current approached, but did not fully achieve, exact symmetry. 

We compare real-time Dst with the orbit-averaged bB N for this storm, as shown in Figure 
11. The top panel shows real-time Dst (the red line) and the orbit-averaged 6 B N (black stars) 
plotted as functions of UT across the entire storm interval. The bottom panel contains the scatter 
plot of orbit-averaged 6 Bn versus Dst. Superposed on the plot are the numerical and graphic (red 
line) results of linear regression analyses performed on the plotted data. The dashed line with a 
slope of unity is provided for reference. For this case, the linear regression slope is near unity 
(0.977) and the correlation coefficient is very high (0.970). It is also clear that the orbit-averaged 
6 Bn data points generally fell below Dst traces. This was most prevalent near storm time 
maximum epochs. It is also reflected in the -9.7 nT intercept obtained through linear regression 
analysis. There are two reasons for the baseline differences. First, the real-time Dst has known 
offsets from final Dst index. Second, such a difference is expected even with final Dst since the 
Dst does not consider stable magnetospheric fields such as the 8 nT from the magnetotail 
currents and a few nT from the quiet time ring current [e.g., Ltihr and Maus, 2010], This example 
along with the others in Le et al. [2011] demonstrate that we can extract a parameter 6 B N whose 
orbit-averaged characteristics mimic those of the provisional Dst index, an important input 
parameter for geomagnetic modeling. 
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5. Concluding Remarks 


The external currents driven by the solar wind-magnetosphere interaction are very 
dynamic and change in various time scales much shorter than those of internal sources. In 
addition, the external currents are primarily ordered by the local time, which is very different 
from the Earth’s internal field. We present recent observations of the Earth’s magnetic field from 
LEO satellites, including polar orbiting ST-5 spacecraft and low-inclination C/NOFS spacecraft, 
and demonstrate that data from multiple spacecraft are required to characterize these external 
currents. Based on our recent observations as well as previous work in literatures, we summarize 
the findings based on these measurements: 

(1) Simultaneous multi-point measurements along a single LEO polar orbit can reveal the 
temporal variability of field-aligned currents in various time scales, measure the motion 
of large scale current sheets, provide opportunities for magnetic gradiometer 
determination of the current density, quantify the closure path of ionospheric Pedersen 
current, and assess the strength of auroral electrojets. 

(2) Measurements from a single equatorial LEO satellite can specify the ring current’s 
temporal evolution, quantify its local time asymmetry, and extract a timely proxy for the 
provisional Dst index at high cadence. 

(3) Field-aligned currents have very complex structures with filamentary currents in various 
scales embedded within large-scale current sheets. Simultaneous measurements with 
longitudinal separations less than ~ 500 km are also required to specify their meso-scale 
variations. 
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(4) Simultaneous monitoring of the dawn-dusk, day-night, and north-south auroral zones are 
also needed to specify the global distribution of field-aligned currents and ionospheric 
currents. This requires placing multiple satellites in polar orbits with large local time 
separations (~ 3 - 6 hours). 

In geomagnetism investigations, great advances have been made since space-based 
magnetic field measurements from dedicated geomagnetism satellites (Magsat, Orsted and 
CHAMP) became available. The upcoming Swarm mission will be the first geomagnetism 
constellation [Friis-Christensen et al., 2006], It contains three satellites, two at lower altitude 
flying side-by-side and one at higher altitude slowly drifting away from the lower-altitude pair 
longitudinally. It will return the first simultaneous geomagnetism measurements at different 
latitudes and longitudes. For post-Swarm geomagnetism satellite missions, it is desirable to have 
a constellation of more than three satellites in a combination of both low and high inclination 
orbits. The constellation would provide simultaneous measurements not only at different 
latitudes and local times, but also with a global coverage. These measurements will result in a 
global specification of the external currents and enable us to separate their magnetic 
contributions from the main field measurements. 

In the post-Swarm era, a desired geomagnetism constellation mission would contain both 
low and high inclination satellites. From the point of view optimal for measuring the external 
currents, the low-inclination satellites would be dedicated to the low-latitude current systems (the 
ring current, the magnetopause current, and the tail current), while high-inclination satellites to 
the high-latitude current systems (the combined field-aligned/Pedersen currents and auroral 
electrojets). Based on our recent observations, we would recommend the following constellation 
configuration: 
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(1) Two or three satellites in the same polar orbit to measure field-aligned currents and 
their temporal variability; 

(2) Additional two or three spacecraft in polar orbits, equally spaced in local time among 
all the polar orbits, to provide the global coverage of magnetospheric and ionospheric 
currents; 

(3) One satellite in a low-inclination orbit to monitor the symmetric and asymmetric parts 
of the ring current. 

Such a constellation is able to provide unprecedented geomagnetism data set with 
simultaneous measurements of current systems at various temporal and spatial scales, 
simultaneous measurements in both in northern and southern polar regions, high accuracy 
and high precision measurements with repeated paths. It allows distinguishing quantitatively 
the external effects from the main internal field at a time scale shorter than an orbit period of 
LEO satellites. 
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(b) 


Figure 1 . Illustration of large-scale electric current systems in the magnetosphere and the 
ionosphere, (a) Magnetospheric current systems; (b) Ionospheric current systems 
(Modified from Le et al. [2010]). 
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Space Technology 5 



Figure 2. Overview of Space Technology 5 (ST-5), NASA’s first three-satellite constellation 


mission in LEO polar orbit (Adapted from Slavin et al. [2008]). 
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Figure 3. Two examples of FAC current density determination using the magnetic gradiometry 
technique. The purple and green traces are from ST-5 pairs 094-224 and 094-155, respectively. 
The black trace shows the traditional single spacecraft FAC current determination applied to 12 s 
averaged data from S/C 094. (Adapted from Slavin et al. [2008]) 
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Spacecraft Footprints in Magnetic Coordinates 
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Figure 4. Overview of field- aligned current observations during the intense magnetic storm on 14 
April 2006. (top) The spacecraft ionospheric footprints for a northern and a southern polar cap 
passes in geomagnetic coordinates. The footprints for the southern pass have been flipped to the 
northern polar cap. (bottom) Overview of ST5 magnetic field variations generated by the field- 
aligned currents. The three components of the magnetic field residual vector (data with the 
internal IGRF model magnetic field removed) are shown in solar magnetic (SM) coordinates. 
The labels for the spacecraft positions (altitudes, magnetic latitudes, and magnetic local times) 
on the bottom are for midspacecraft 094 (black) only. (Adapted from Le et al. [2009]) 
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FACs during Northern Polar Cap Pass 
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Figure 5. The time-shifted FAC magnetic field variations displayed in the MVA coordinate 
system for the northern polar cap pass during the intense magnetic storm on 14 April 2006. The 
arrows indicate the general direction of the large-scale currents with a downward arrow for 
currents flowing into the ionosphere and an upward arrow for currents flowing out of the 
ionosphere. (Adapted from Le et al. [2009]) 
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Figure 6. The FAC current setup and geometry for simple calculations of the magnetic field 
signatures, (a) Three pairs of (left) balanced infi- nite current sheets to model the two pairs of 
(right) unbalanced current sheets, (b— d) Current (left) density distributions and their (right) 
magnetic field signatures. (Adapted from Le et al. [2010]) 
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Figure 7. Two examples of ST-5 polar cap crossings. The bottom panels show the magnetic field 
signatures of these currents. The top panels show the least-square fit current density distributions 
along the orbit track. In the bottom panel, the red traces are the best fit magnetic fields, and the 
black traces are the observed magnetic fields in the cross-track direction along the Y axis. 
(Adapted from Le et al. [2010]) 
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Figure 8. Statistic results of the R2 current intensity versus the R1 current intensity. The solid 
line has a slope, which is the average of the R2 intensity to R1 intensity ratio. The dashed line 
has a slope of 1 . (Adapted from Le et al. [2010]) 
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Figure 9. Overview of the July 22, 2009 storm containing 5 days of hourly averaged OMNI data 
from July 21 to 25, 2009. Shown from top to bottom are the interplanetary magnetic field (IMF) 
Bz component, the solar wind density, velocity, dynamic pressure, and the Kyoto provisional Dst 
index. (Adapted from Le et al. [2011]) 
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Figure 10. Local time variations of <5Bn at various stages of the storm evolution, corresponding 
to the times marked by the dashed lines in Figure 10. (Adapted from Le et al. [2011]) 
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Figure 11. Comparison between the Kyoto Dst index and the provisional Dst index determined 
from the C/NOFS magnetic field data for the July 22, 2009 magnetic storm. (Adapted from Le et 
al. [2011]) 
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